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Abstract 

We examine the kinetics of a snowball that is gaining mass while is 
rolling downhill. This dynamical system combines rotational effects with 
effects involving the variation of mass. In order to understand the con- 
sequences of both effects we compare its behavior with the one of some 
objects in which such effects are absent, so we compare the snowball with 
a ball with no mass variation and with a skier with no mass variation nor 
rotational effects. Environmental conditions are also included. We con- 
clude that the comparative velocity of the snowball respect to the other 
objects is particularly sensitive to the hill profile and also depend on some 
retardation factors such as the friction, the drag force, the rotation, and 
the increment of mass (inertia). We emphasize that the increase of inertia 
could surprisingly diminish the retardation effect owing to the drag force. 
Additionally, when an exponential trajectory is assumed, the maximum 
velocity for the snowball can be reached at an intermediate step of the 
trip. 



1 Introduction 

The snowball is a dynamical object that gains mass while is rolling downhill. It 
is a particularly interesting problem since it permits to combine several concepts 
of mechanics such as the traslational and rotational dynamics, the rigid body, 
mass variable systems, normal and tangential forces, radius of curvature for a 
trajectory etc ^G|- Modeling this problem implies many input parameters and 
the use of numerical methods. Additionally, an ansatz should be made about 
the way in which the mass (or volume) is growing with time. Environmental 
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conditions are also considered utilizing well known assumptions about friction 
and drag forces. The dynamical behavior of the snowball will be studied in the 
velocity vs time, and velocity vs length —planes. 

Moreover, comparison with other dynamical objects could clarify many as- 
pects of the very complex behavior of the snowball. Therefore, we will develop 
our analysis by comparing the snowball (SB) motion with the one obtained 
from a skier sliding without friction (SNF), a skier sliding with friction (SF) 
and a ball with constant mass and volume (B). 

In section|21 we discuss the basic assumptions and write out the equations of 
motion for the snowball. In section |21 the comparison between the four dynam- 
ical objects mentioned above is performed in the asymptotic regime. Section 
^describes some proves of consistency to test our results. Section |S] provides 
a complete analysis of the comparative kinetics of the four dynamical objects 
based on some environmental retardation factors. Section El is regarded for the 
conclusions. 



2 The problem of the snowball rolling downhill 

2.1 Basic assumptions 

The complete analysis of a snowball requires many input parameters. The 
problem does not have any analytical solution, so it has to be solved numerically 
We have listed the numerical values assumed in this paper in table [2 on page 
00 Besides, we make the following assumptions in order to get an equation of 
motion 

1. The snowball is always spherical meanwhile is growing and acquiring mass. 
Its density is constant as well. 

2. It is supposed that there is only one point of contact between the snowball 
and the ground, and that the snowball rolls without slipping throughout 
the whole motion. Wc shall call it the No Slipping Condition (NSC). 

3. In order to accomplish the previous condition, the frictional static force 
that produces the rotation of the snowball (Fr s ), has to hold the condition 
Fr s < where N is the normal (contact) force to the surface, and /i s is 
the coefficient of static friction. We assume that /i s is independent on the 
load (Amontons' law), this statement applies for not very large loads 0]. 

4. The drag force owing to the wind is assumed to be of the form 

Fv = -e*M* (i) 

where is the air density, Cd is the air drag coefficient, and A is the 
snowball's projected frontal area i.e. A = nr 2 . 

We assume the air drag coefficient Cd to be constant, since this assumption 
has given reasonable results in similar problems On the other hand, 
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it has been established that the force F v could be a linear or quadratic 
function of the speed depending on the Reynolds number (Re) For 
Re > 1 (which is our case) a quadratic dependence fits well with experi- 
mental issues jH]- 

5. The mass of the snowball increases but finally reaches an asymptotic value. 
Furthermore, a specific function of the mass (or volume) in terms of time 
must be supposed. In our case we assume the following functional form 



where Mq is the initial mass of the snowball and clearly the final mass is 
Mo + Kq. 

6. A hill profile must be chosen in order to make the simulation, like in other 
similar problems [3] , it has an important effect on exit speed. Specifically, 
we have chosen an inclined plane and an exponential trajectory. 

2.2 Equations of motion 



Figure 1: A snowball rolling downward on a wedge, the X — axis is parallel to 
the wedge surface and the Y — axis is perpendicular. 

To simplify the problem we start assuming the snowball rolling downhill 
on a wedge whose angle is 8 (see Fig. In the time t the snowball has a 
mass M and in the time t + dt its mass is M + dM, let us consider a system 
consisting of the original snowball of mass M plus the little piece of snow with 




(2) 




4 



R. A. Diaz, D. L. Gonzalez, F. Marin, and R. Martinez 



mass dM. At the time t, the momentum of such system is P (i) = Mv (bold 
letters represent vectors) since the little piece of snow is still on the ground and 
at rest. At the time t + dt the ball has absorbed the piece of snow completely, so 
the momentum is given by P (t + dt) = (M + dM) (v+dv) then the momentum 
change is dP = Mdv + vdM (where we have neglected the differential of second 
order) and the total force will be 

„ dP rfv dM 

F = — =M— +v— 3) 

dt dt dt w 

where v corresponds to the velocity of the center of mass respect to the ground. 

Further, the system is rotating too, we also suppose that such rotation is 
always made around an axis which passes through the center of mass (perpen- 
dicular to the sheet in Fig. In this case the rotation is around a principal 
axis, hence the equation of motion for the angular momentum is given by 

L c = Ic~^c (4) 

where the subscript C refers to the center of mass coordinates. Ic denotes the 
moment of inertia of the snowball measured from an axis passing through the 
center of mass, and ~uj c refers to the angular velocity. According to figure 
Lc is directed inside the sheet, and the torque will be 

dL c _> . . 

-dT =Tc (5) 

Wc should remember that this equation is valid even if the center of mass is 
not an inertial frame PP, which is clearly our case. To calculate dLc we make an 
analogous procedure as in the case of dP, and the equation JSJ is transformed 
into 

d~uj c . -> dl c -> 

Ic — + UJC ^r= TC - (6) 

where is the total torque measured from the center of mass. For the sake 
of simplicity, we will omit the subscript C from now on. 

On the other hand, the external forces and torques exerted on the system 
are similar to the ones in the simple problem of a ball on a wedge PP 

F = W + N + F Rs + F a ; ~t = r x F Rs (7) 

where W is the weight (which acts as the driving force), N the normal force, 
Fjis the statical friction force, and F a is any applied force which does not 
produce torque. If we use as a coordinate system the one indicated in figure ^ 
(with the Z — axis perpendicular to the sheet) wc can convert these equations 
into scalar ones; then using Eqs. © and Q we get by supposing that F a is 
parallel to the X — axis. 
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The first equation decouples from the others whenever the NSC is maintained; 
so we will forget it by now. We should be careful in utilizing the NSC since the 
radius is not constant. The correct NSC in this case consists of the relation 
ds = r (dip) where ds is the differential length traveled by the center of mass of 
the snowball (in a certain differential time interval dt), d<p is the angle swept 
along the differential path, and r is the snowball radius in the time interval 
[t,t + dt]. Using the correct NSC, we get 



ds 
dt 



dt 



(9) 



and taking into account that the radius depends explicitly on the time, we obtain 

(10) 



a 1 dr 

a = — lj 

r r dt 



where uj is the angular velocity, a is the angular acceleration, and v, a are the 
traslational velocity and acceleration respectively 1 . 

It is convenient to write everything in terms of the displacement s, taking 
into account the following relations 



ds d s 

v =dt:> a =d£ ; 

replacing pfT?HfTT|l into ©, we find 

~1 d 2 s 



dt 



rF Ra = I 



dt 2 



1 dr ds 
^~dt~dt 



1 ds 
rdt 

1 ds 



dl 
~dt 



(11) 



(12) 



Now we use the moment of inertia for the sphere / = (2/5) Mr 2 and the fact 
the the Mass "M" and the radius "r" are variable, then 



dl 2 / , dM _ dr 



2Mr- , (13) 



dt 5 V dt dt 

4 , dM 9 dr 
M = -irpr 3 ; — =^(yr 2 — (14) 



1 Observe that we can start from the traditional NSC with v = cur. Nevertheless, the other 
traditional NSC a = ra, is not valid anymore. 
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where the snowball density p has been taken constant. Additionally, we assume 
the applied force F a to be the drag force in Eq. (JTJ. With all these ingredients 
and replacing Eq. (|12[) into equations JSJ, they become 



as 



l5 PA C d 1 
56p r 



23 1 dr ds 
Tr~dldt, 
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(15) 



Finally, to complete the implementation of Eqs. I|15f) . we should propose an 
specific behavior for r (i) or M (t) ; we shall assume that M (t) behaves like in 
Eq. (J2J). On the other hand, as a proof of consistency for Eqs. (|15fl . we see that 
if we take r = constant, and Cd = 0, we find that 



d 2 s 5 

dfi = 7 9Smd 

F Rs = ^M g sinO (16) 

which coincides with the typical result obtained in all common texts for the ball 
rolling downward on a wedge with the NSC [J. 



2.3 A snowball downhill on an arbitrary trajectory 




Figure 2: A snowball rolling downward on an exponential trajectory. To find the 
local value of the angle 9 we define the X, Y — axis as indicated in the figure. 
We see that for a sufficiently large value of x we get 9^0. 



In this case the acceleration written above converts into the tangential one, 
and the equation (JHJ) for the normal force becomes 

v 2 

N - Mgcos9 = M— (17) 
R 
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where R is the radius of curvature. In order to solve Eq. i|17|). it is convenient to 
use the coordinate axis plotted in figure |2 In cartesian coordinates, the radius 
of curvature is given by 



R(x) = 



3/2 



y" 
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( llL 

dx 



(18) 



Moreover, the angle 9 is not constant any more, and according to the figure [21 
we see that 



sm ( 



—dy 
ds 

dx 
ds = 



—dy 
(dxf + (dy) 2 



1 



i + (y'Y 



i + (y'Y 



(19) 



where the minus sign in the differential dy is due to the decrease of the coordinate 

y- 

So the problem of the snowball rolling on an arbitrary trajectory can be 
solved by replacing l|19|) into the first of Eqs. 1|15[) . and making an assumption 
like (J2J). Additionally, Ecr. (|17|) provides the solutions for the normal force by 
considering Eqs. i|18|) . (|19|) and the solution for the velocity obtained from (|15f) . 
Notwithstanding, the first of Eqs. (|15|) docs not depend on the normal force 
whenever the NSC is maintained. So, we shall ignore it henceforth. 



3 Comparison of the snowball with other dy- 
namical objects 

May be the clearest way to understand the dynamical behavior of the snowball, 
is by comparing it with the dynamics of other simpler objects. In our case we 
shall compare the dynamics of four dynamical objects 

1. A skier sliding with no friction (SNF). 

2. The same skier but sliding with friction (SF). 

3. A ball with constant mass and volume (B) 

4. A snowball with variable mass and volume (SB) 

Such comparison will be performed in the v — t and v — x planes. The 
behavior of the first two objects were reproduced from the article by Catalfamo 
0, and the equations for the ball were obtained from the ones of the snowball 
by setting r — > constant. In making the comparison we use the input parameters 
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M = K = 85 


p = 917 


PA = 0.9 


Us = 0.03 


C d = 0.3 


= (4.76) u 


A s = 0.6 


/3 = 0.07 


g = 9.8 


h = 25, a = 0.035 


Ho = 0.03 


v = 



Table 1: Input parameters to solve the equation of motion for the SNF, SF, 
SB, and the B. All measurements are in the MKS system of units. Mo is the 
initial mass of all objects, Ko defines the increment of mass for the SB see eq.(J2J), 
p and pa are the snow and air densities respectively, p, s the statical coefficient 
of friction between the B and the ground (and also between the SB and the 
ground). Cd is the air drag coefficient, 9 the angle of the wedge, As is the skier 
frontal area, /3 is a parameter that defines the rapidity of increase of mass in the 
SB sec Eq. (J3J). g is the gravity acceleration, /io is the dynamical coefficient 
of friction between the SF and the ground, and vo is the initial velocity of the 
four objects. Further, in the exponential trajectory y = he~ ax , where h is the 
height of the hill. 

of tabled] most of the parameters in this table were extracted from P] 2 . As it 
can be seen, we assume the same mass Mo for the skier and the ball, and this 
is also the initial mass of the snowball ending with a mass of 2Mo- 

3.1 Asymptotic limits 

Before solving all these problems, we shall study the asymptotic limits of the four 
objects in the inclined plane and the exponential trajectory with and without 
drag force. The asymptotic regime provides useful information and can be used 
to analyze the consistency of the numerical solutions. These limits depend on 
the drag force, the trajectory and the object itself. 

3.1.1 Inclined plane with no drag force 

For each object we obtain the following limits 

• For the SF its velocity is easily found 

v — vo + gt (sin 9 — [Ib cos 9) (20) 

so there is no finite value for the velocity and its behavior is linear respect 
to time. The SNF asymptotic limit is obtained just taking /jld — ► 0. 

• For the SB from Ea. (|15[l and assuming that the radius reaches an asymp- 
totic limit i.e. dr/dt — > when t — > oo we get 

5 

v (t — > oo) = = Vq + —gtsin9 (21) 

2 However, we have changed mildly the parameters that describe the profile of the expo- 
nential trajectory, in order to keep the NSC throughout. 
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getting a linear behavior respect to time. The same behavior is exhibited 
by the B when t — > oo. 

Observe that Voo m the SB is independent on the mass and equal to the 
value for the ball B, it is reasonable owing to the asymptotic behavior 
assumed for the SB, Eq. 

3.1.2 Inclined plane with drag force 

• For the SF its equation of motion is easily obtained 

dv PACdA 2 

— = -nog cos9 v + gs\nO 

and in the asymptotic limit dvoo/dt — > 0, so we get that 

t& = — ^^-(sinfl-^cosfl) (22) 
PA^d A 

Now, by setting n, D — > 0, we get for the SNF 

• For the SB from Eq. (fT5|) and setting d 2 s/dt 2 — ► 0, dr/dt — ► 0, we obtain 

Voo 

u oo = Tf, pr?- oo 0sin0 = — —-—sinS (23) 

where r M = r (t — > oo) . The second term in <|23[) is obtained from Eq.J2J by 
taking the asymptotic limit when t — > cxd, in this case is given by 

7^- (M + K ) 

In the case of the B, we obtain the expression (|23|) but r^, M 00l are 
constant in time and equal to its initial values To, Mq, Aq. 

3.1.3 Exponential trajectory with no drag force 

In this case it is easier to examine the asymptotic limits, since when the objects 
have traveled a long path, 6 — » and the object is reduced to run over a 
horizontal plane, see figure |3 Therefore the limits are 

• For the SF, = 0. 

• For the SNF, it is found easily by energy conservation 

vlo = 2.9^ + t>0 



where h is the height of the hill. 
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• For the B we can find Voo by taking into account that energy is conserved 
because friction does not dissipate energy when the NSC is held pQ. By 
using energy conservation we obtain 3 

For the snowball the limit is not easy to obtain because energy is not con- 
served and Eq. (|15f) does not provide any useful information 4 . However, 
according to the figures [5] this limit is lower than the one of the B owing 
to the increment of inertia. 

3.1.4 Exponential trajectory with drag force 

In this case all velocities vanish for sufficiently long time. 



4 Proves of consistency 

The equations of motion for each object where solved by a fourth order Rungc 
Kutta integrator[7]. To verify the correct implementation of the program, wc 
reproduce the results obtained by and solve analitically the problem of the 
ball of constant radius in the inclined plane without drag force, the results were 
compared with the numerical solution. Additionally, all the asymptotic values 
discussed above were obtained consistently. 

Finally, the reader could have noticed that one of our main assumptions 
was the NSC. However, this condition can only be valid if we ensure that the 
statical friction force does not exceed the value /j, s N all over the trajectory in 
each case. Otherwise, the snowball would start slipping preventing us of using 
the fundamental relations and (|10|l . Additionally, if the snowball started 
sliding, the frictional force would become dynamics i.e. Fr d = /idN producing 
a coupling among the three Eqs. ©, remember that we have assumed that the 
first one was decoupled. As a first approach, we study the validity of the NSC 
in the asymptotic regime by calculating Fr s (t — > oo) and fJ, s N (t — > oo). 

For the inclined plane, these limits can be estimated by using Eq. @, the 
first of Eqs. ©, and the second of Eqs. lfT3|) 

/j. s N (t — > oo) — > fisMocgcosO 

2 fd 2 s\ , x 

F Rs (t^^) - 5 M »(^fJ ■ ( 24 ) 

it is easy to verify that the NSC is valid in the asymptotic limit for the wedge, 
since in the case of the presence of a drag force we get Fr s (t — > oo) = 0. Addi- 
tionally, in the case of absence of the drag force it is found that Fr s (t — > oo) = 

3 Since what really matters for energy conservation is the height of the center of mass, there 
is a tiny difference that can be neglected if the radius of the ball is much smaller that the 
height of the hill. 

4 By using 9 — > 0, and d 2 s/dt 2 — > 0, the second of Eqs. 1151 becomes trivial. 
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(2/7) (Mo + K ) gsin8, and the condition Fr s (t — > oo) < fj, s N (t — > oo) is ac- 
complished by our input parameters (see table 0). 

For the exponential trajectory the analysis is even simpler, since the path for 
large times becomes a horizontal straight line, the asymptotic limits for Fr s and 
fi s N are the same as in Eqs. <|24ll but with 9 — > 0; and the NSC condition is 
held when t — > oo. 

However, the NSC in the asymptotic regime does not guarantee that it is 
held throughout the path. For example, in the case of the exponential trajectory, 
the maximum slope of the profile is found at the beginning of the trajectory, it 
was because of this fact that we changed mildly the profile parameters defined 
in Ref. |3J. Consequently, by using the first of Eqs. © as well as the Eqs. (|15fl. 
l(T7|) . l(TH|) . and (|T§|) ; we solved numerically for Fr s vs time and length and for 
/i s N vs time and length; utilizing the numerical input values of tabled We then 
checked that Fr s < fi s N throughout the time or length interval considered in 
each case. 

5 Analysis 

In order to make a clear comparison, we take the initial mass of all the objects 
to be equal, and all initial velocities are zero. In figures ISTol we use the following 
conventions: The dashed dotted line corresponds to the SNF, the solid line 
represents the SF, the dashed line refers to the B, and finally the dotted line 
corresponds to the SB. In both, ball and snowball we only consider statical 
friction and neglect possible dynamic frictional effects due to sliding, because we 
have guaranteed that the NSC is valid throughout the trajectory, as explained 
in section 01 



In figure [3] we plot d vs t and v vs x for constant slope of the wedge 
without drag force. Of course, all graphics in the v — t plane are straigth lines 
except the one for the snowball. We can see that the SNF is the fastest object 
as expected, since no retardation factors are acting on it, next we have the B 
which posseses the rotation as a retardation factor. Additionally, the SB line 
is always below the B line because in the former, two retardation factors are 
present: the rotation and the increase of inertia. However, for sufficiently long 
time (or length) the increase of inertia vanishes (according to our assumptions) 
so that the velocities of both B and SB coincide, in agreement with the analysis 
made in Sec. 13.1.11 We checked that condition, though it does not appear in 
the Fig. [3J because of the short time and length interval displayed. 

The line corresponding to the SF is below the line corresponding to the 
SNF as it must be, however the relation between the SF and SB lines is 
particularly interesting and deserves more attention. At the beginning the SB 
is slightly slower than the SF, but for sufficiently long time, the SB becomes 
faster. It can be explained in the following way, at the beginning the SB has 
two retardation factors: the rotation and the increase of inertia, while the SF 
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Figure 3: Plots in the v — t plane (left) and the v — x plane (right) when the 
objects travel in a wedge of constant slope with no drag force. 



only has one retardation factor: the sliding friction. On the other hand, for 
sufficiently long time the increment of inertia becomes negligible in the SB, and 
only the rotation acts as a retardation factor, consequently the SB behaves like 
the B as shown in section 13 Therefore, the combination of the two retardation 
factors at the beginning makes the SB slower but when the increase of inertia 
is small enough, the SB becomes faster than the SF. Nevertheless, we should 
point out that this behavior depend on the value of [Id, if it were large enough 
the line for the SF would lie below the lines for B and SB at all times, in 
contrast if it were small enough the SF line would lie above the B and SB lines. 
Notwithstanding, the rapidity of the SF must be smaller than the SNF speed 
at any time and for any value of fj,n- 

According to this analysis, when the objects travel in a wedge with no drag 
force, the pattern of velocities in descendent order for any set of the input 
parameters (as long as the initial masses and velocities are the same) is the 
following: the SNF, the B and the SB. The comparative velocity of the SF 
depend on the input parameters but it must be always slower than the SNF. 
As a proof of consistency, it can be checked that the asymptotic limits in Eqs. 
iffiH) , obey this pattern. 
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Figure 4: Plots in the v — t plane (left) and the v — x plane (right) when the 
objects run over a wedge of constant slope with drag force. 



Figure |4| correspond to a wedge with constant slope including drag force. 
In this case the comparative behavior among the four elements is not as simple 
as in figure because in this case the lines cross out each other. However, the 
line describing the SF is always below the line describing the SNF as it should 
be. This more complex behaviour owes to the frontal area dependence of the 
drag force. For instance, we can realize that at short times the comparative 
behavior is very similar to the one in figure [3] since the drag force has not still 
acted significantly. All these elements get an asymptotic limit as we described 
in section 13 We see that the largest asymptotic limit correspond to the SB, 
in opposition to the case of figure with no drag force, in which the snowball 
was one of the slowest objects; the clue to explain this fact recides in the frontal 
area dependence of the drag force. From Eqs. (|22l I23|) we can verify that for 
all these objects the terminal velocity behaves as v 2 oc M/A, this quotient is 
larger for the B than for the SNF and the SF in our case, then the asymptotic 
velocity vb is larger than vsnf and vsf, for both the skier and the ball this 
ratio is a constant. In contrast, since in the snowball the mass grows cubically 
with the radius while the area grows quadratically, its velocity behaves such 
that Vg B oc r (t) . Therefore, for sufficiently long times, its velocity grows with 
the radius of the SB, getting a higher terminal velocity (of course it depends on 
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Figure 5: Plots in the v — t plane (left) and the v — x plane (right) when the 
objects travel on an exponential trajectory with no drag force. 



the asymptotic value of r (t)). Observe that if we had assumed a non asymptotic 
behavior of r (t) in Ijl5|l we would have not obtained any finite terminal velocity 
for the snowball even in the presence of a drag force. Furthermore, we see that 
the terminal velocity for the SB is reached in a longer time than the others, it 
is also because of the slow growth of r{t). 

Figure [5] describes the four elements traveling in an exponential hill with 
no drag force. Two features deserve special attention: (1) the terminal velocity 
is achieved in a very short time specially in the cases of the SNF and the B, 
these limits coincides with the ones obtained in section 01 (2) For the SB and 
the SF there is a local maximum velocity at a rather short time, the diminution 
in the velocity since then on, owes to the decreasing in the slope of the path 
in both cases, the increment of inertia in the case of the SB, and the friction 
in the SF. Such local maximal velocity cannot be exhibited by the SNF and 
the B because conservation of energy applies for them, and as they are always 
descending their velocities are always increasing, though for long times they 
are practically at the same height henceforth, getting the terminal velocity. In 
particular, we see that the terminal velocity of the SF is zero as it was shown 
in Sec. ETT^l 
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Figure 6: Plot in the v — t plane (left) and the v — x plane (right) when the 
objects run over an exponential trajectory with drag force. All velocities vanish 
for sufficiently long time (or length) as expected, despite it is not shown in the 
interval displayed. 

In figure |6] the elements travel in an exponential hill with drag force. In 
this case, the conservation of energy does not hold for any of the objects, conse- 
quently maximum velocities in intermediate steps of the trajectory are allowed 
for all of them. All terminal velocities are zero as expected. Because of the 
same arguments discussed above, the line of the SF is below to the one of the 
SNF. However, any other pattern depend on the numerical factors utilized. 

A final comment is in order, we can realize that though the solution of the 
kinetics of the snowball depends on the ansatz made about the mass growth, 
the bulk of our results and analysis only depend on the fact that the snowball 
mass reaches a finite asymptotic value. So that the discussion is quite general, 
especially in the asymptotic regime. 

6 Conclusions 

We have described the behavior of a snowball acquiring mass while rolling down- 
hill, taking into account the enviromental conditions. The dynamics of the 
snowball is very complex because it is a rotating object and at the same time 
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its mass and moment of inertia are variables. In order to visualize better the 
effects due to the rotation and mass variation, we compare its motion with the 
kinetics of two objects in which the rotation and mass variational effects are 
absent (the Skier with Friction and the Skier with No Friction) , and with one 
object in which the rotation is present but no the mass variation (the Ball of 
constant mass and radius). 

The comparative behavior of these objects depend on the trajectory but 
also on some retardation factors: the friction, the drag force, the increase of 
mass (inertia), and the rotational effects. It worths to remark that despite the 
increment of inertia is a retardation factor in some circumstances, it could sur- 
prisingly diminish the retardation effect due to the drag force. In addition, some 
local maxima of the velocities for each object appears in an exponential trajec- 
tory, showing that the maximum velocity might be achieved at an intermediate 
step of the path. 

Finally, we point out that despite the complete solution of the snowball de- 
pends on an ansatz about the way in which its mass grows; its comparative 
dynamics respect to the other objects in the asymptotic regime is basically in- 
dependent of the details of the growth, and only depend on the assumption that 
the mass reaches an asymptotic value, a very reasonable supposition. Therefore, 
we consider that our analysis is not very model dependent at least in the regime 
of large times or lengths. In addition, these asymptotic limits serves also to 
show the consistency of our results. 
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